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THE  DYNAMIC  VISCOPLASTIC  EXPANSION 


OF  A  CYLINDRICAL  TUBE 


by 


Edmund  J.  Appleby 


ABSTRACT 

The  viscoplastic  flow  of  a  long  thick-walled  tube  is  in¬ 
vestigated.  The  tube  is  subjected  to  internal  pressure  and  has  its 
ends  restrained  from  motion  in  die  axial  direction.  The  material 
of  the  tube  is  rigid-viscoplastic  and  incompressible.  The  pressure 
required  to  produce  a  specified  expansion  of  die  tube  is  calculated 
for  two  examines.  In  tbs  former  die  effect  of  different  viscosity 
coefficients  is  observed*  In  die  latter  example  a  comparison  is 
made  of  the  effects  of  perfect  plasticity,  viscosity  and  inertia. 


1.  INTRODUCTION 


A  study  is  made  of  the  mechanical  behavior  of  a  long  hollow  circular 
cylinder*  with  its  end  restrained  from  motion  in  the  axial  direction*  when 
it  is  subjected  to  internal  pressure.  The  basic  problem  of  plane  strain 
with  rotational  symmetry*  being  one  of  those  which  most  readily  yield  to 
treatment  in  plasticity*  is  a  standard  problem  in  this  field*  and  has  been 
studied  by  investigators  in  a  variety  of  different  ways  for  the  ideally 
plastic  and  the  elastic  plastic  material.  Reviews  of  such  solutions  have 
been  given  by  Hill  [  l]  *  and  Prager  and  Hodge  [2].  The  present  paper 
considers  an  ideal  material*  a  viscoplastic  Bingham  solid  [  3]  *  which 
is  undeformable  until  the  stresses  reach  their  yielding  values*  and  then 
under  stresses  which  are  in  excess  of  their  yielding  values,  has  strain 
velocities  dependent  on  this  stress  excess  or  overstress.  The  material 
does  not  exhibit  work  hardening,  and  as  it  is  known  that  in  the  fully 
plastic  state  volume  changes  are  negligible*  it  is  assumed  to  be  incom¬ 
pressible. 

The  solution  is  obtained  by  using  the  viscoplastic  constitutive 
equations  due  to  Hohenemser  and  Prager  [4]  .  The  analysis  is  valid  for 
a  general  plastic  yield  condition.  An  approach  by  means  of  a  linearized 
theory  of  visco plasticity*  in  which  the  flow  is  specified  by  Prager' s  constitu¬ 
tive  equation  [5]  ,  is  equivalent. 

A  general  solution  is  formulated,  but  in  order  to  simplify  the  numeri¬ 
cal  calculations  a  specific  expansion  is  imposed  on  the  tube  in  which  the 
interior  boundary  has  a  uniform  radial  acceleration.  The  pressure-time 
variations  required  to  maintain  this  flow  for  different  viscosity  coefficients 
are  compared.  A  different  expansion  is  then  imposed,  in  which  the  inner 
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radius  of  the  tube  expands  to  one  and  one  half  its  initial  size.  beginning 
and  ending  with  zero  velocity.  Comparison  is  made  of  the  effects  of 
perfect  plasticity,  viscosity  and  inertia.  If  the  expansion  takes  place 
slowly,  the  effect  of  inertia  on  the  required  pressure  is  negligible, 
and  the  effects  of  perfect  plasticity  and  viscosity  are  comparable.  On 
the  other  hand  if  the  expansion  occurs  very  quickly  the  effect  of  inertia 
becomes  comparable  with  that  of  viscosity,  and  the  perfectly  plastic 
contribution  to  the  pressure  is  negligible. 


2.  BASIC  EQUATIONS 


Let  the  space  variables  be  a  system  of  cylindrical  co-ordinates 

« 

r,  9  ,  z  in  which  the  z-axis  coincides  with  the  axis  of  symmetry(  then 
the  tube  is  bounded  by  the  cylinders  r  =  a  and  r  a  b,  where  a  <  b. 

As  the  pressure  inside  the  tube  increases)  the  stresses  in  the  material 
nearest  the  interior  boundary  will  be  the  first  to  reach  the  yield  limit. 

With  further  increase  of  pressure  the  plastic  region  will  extend  until 
its  outer  boundary  coincides  with  the  tube's  outer  surface.  Until  this 
state  is  reached)  the  flow  of  the  plastic  innermost  region  of  the  tube 
is  restricted  by  a  surrounding  rigid  region  and  by  the  conditions  of  axi¬ 
ally  symmetric  plane  strain  on  the  incompressible  material.  The  whole 
tube  therefore  remains  rigid)  and  hence  the  stresses  in  the  innermost 
plastic  region  reach  but  do  not  exceed  the  yield  limit.  When  the  whole  of 
the  tube  becomes  plastic(  the  material  is  about  to  flow  in  an  unrestricted 
manner,  since  any  further  increase  in  the  internal  pressure  will  then 
produce  overstress  in  the  material  of  the  tube.  The  time  t  is  measured 
from  this  instant,  and  the  values  of  a  and  b  for  t  ss  0  are  denoted  by 

ag  and  bg. 

It  is  assumed  that  the  tube  is  sufficiently  long  to  make  the  stresses 
and  strains  independent  of  the  axial  co-ordinate,  and  that  at  any  instant  of 
the  flow  process  each  particle  of  the  tube  is  moving  radially  outward 
with  a  velocity  u  which  depends  only  on  the  radial  co-ordinate  r.  The 
velocity  components  at  a  point  distance  r  from  the  axis  at  time  t  are 
then 

ur  =  u(r,  t),  ue  =  0  ,  uz  =  0.  (1) 
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The  strain  rates  are 
V  Ju 

*r"ff7  * 


W  •  V0’ 


%z  =  0' 


Vzr  =  0  *  %0=0 


(2) 


where  the  dot  denotes  di£ferentiation  with  respect  to  the  time. 
Since  the  material  is  assumed  incompressible 


8  u 

5~T 


+  -  a  0. 


(3) 


This  differential  equation  has  the  solution 


u  = 


i(t) 

r 


(4) 


e 

in  which  *{t)  is  an  arbitrary  function  of  the  time.  Expressing  the 

velocity  u  of  a  particle(  distance  r  from  the  axis  at  time  t  as 

dr/dt,  and  integrating  (4),  we  find  that  a  particle  initially  at  distance 

I — 2 - 

rg  from  the  axis  is  at  distance  vrg  +  2*(t)  after  a  time  t. 

The  radial  and  circumferential  strain  rates  can  now  be  written 


r  r 


By  the  rotational  symmetry  of  the  flow  fieldt  the  shearing 
stresses  with  respect  to  the  cylindrical  co-ordinates  are  zero,  so 
the  only  equation  of  motion  which  is  not  identically  satisfied  is 


8c 

r 
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=  D 


du 
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(6) 


where  D  is  the  density  of  the  material.  When  u  is  replaced  by  the 
expression  in  (4),  the  above  equation  can  be  written 
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(7) 


The  boundary  conditions  throughout  the  yielding  process  are 
<rr  =  -  p  at  r  =  a, 
ar  =  o  at  r  =  b; 


(8) 


it  being  assumed  that  there  is  no  pressure  on  the  external  boundary* 


The  analysis  so  far  is  independent  of  the  constitutive  equations, 
and  is  the  same  for  all  materials. 

For  any  incompressible  isotropic  material  in  plane  strain  Geiringer 
[  6]  has  shown  that  a  general  yield  condition  can  be  expressed  as  a 
function  of  the  single  variable  -  a a }  and  a2  being  the  principle 
stresses  in  the  plane  z  =  const.  For  the  rotationally  symmetric  pro¬ 
blem  the  principle  stresses  are  a  and  a0  ,  so  that  the  yield  function 
F  can  always  be  written  in  the  form 

F  =  la0  -  orr  1  -  21c  ,  (9) 

where  k  is  the  yield  stress  in  shear*  In  particular  if 

an  assumption  which  may  be  verified  a  posterior,  (9)  can  be  replaced 

by 


F=a0  -  ar  -  2k  (10) 

The  viscoplastic  flow  rule  of  Hohenemser  and  Prager  can  be  written  in 
the  general  form 


X  c , .  =  <F> 


8  F 

sir. 


ij 


(ID 


where  the  notation  is  defined  by: 
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<F>  «  F  if  F  *  0, 

<F>  «  0  if  F  <  0, 

X  is  the  viscosity  coefficient  of  the  material,  and  F  is  the  yield 
function.  For  the  yield  function  ( 10),  therefore,  the  only  nonvanish¬ 
ing  strain  rate  components  for  t  >  0  are  given  by 

X«r  =  -  (<rff  -<rr  -  2k)  ,  X^  »  <r0  -<ry  -  2k.  (12) 

It  is  interesting  to  note  that  since  ( 10)  is  a  linear  function,  ( 12) 
may  also  be  obtained  as  a  special  case  of  the  piecewise-linear  visco¬ 
plastic  flow  law  proposed  by  Prager  [  5]  . 


3.  SOLUTION 


Equations  (12)  can  alternatively  be  written 


a  <ra  -  <r  -  2k. 
0  r 


(13) 


The  substitution  of  equation  (13)  into  the  equation  of  motion  (7)  gives 


(14) 


Integrating  this  equation  with  respect  to  r,  and  employing  the  bound¬ 
ary  condition.  O’  =0  when  r  =  b,  to  evaluate  the  arbitrary  function 
of  time  furnishes  the  stress  components: 


(15) 


Since  a  and  b  can  be  expressed  in  terms  of  their  initial  values 
a.Q  and  bg  and  the  function  *  (t),  equations(  15)  furnish  the  stress 
distribution  at  any  instant  in  terms  of  ♦  and  its  first  and  second 
derivatives.  The  functions  *  and  p  must  satisfy  the  boundary  con¬ 
dition  »  -  p  on  r  ■  a.  Hence,  by  using  the  expression 

r  =  ^rQ2  +  2*, 
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*Dlog 


»20  +  2* 

»20  *  2* 


+ 


0  "  b  0^  +  *2P^b  O'  a  V 

(a2Q  +  2*)  (bZQ  +  2*)  (a20  +  2*)  (b2Q+  2*) 


2p  -  2k  log 


+  2* 
+  2* 


(16) 


the  conditions  at  t  =  0  being  ♦  =  0  and  *  «  0,  that  is 
is  initially  and  the  initial  radial  velocity  is  zero. 


nTT" 
r  0 


+  2* 


If  the  effects  of  inertia  and  viscosity  had  been  neglected* 
the  left  hand  side  of  equation  (16)  would  be  zero*  and  hence  p  =  k  log 
[  1  +  (b2^  -  a2g)/a2]  ,  which  is  the  internal  pressure  required  in  the 
case  of  an  ideally  plastic  tube  to  maintain  it  in  a  state  of  unrestricted 
flow  (cf.  [2]  ,  page  1 18).  If,  on  the  other  hand,  the  inertia  term  alone 
is  neglected  and  p  is  regarded  as  a  constant  internal  pressure, 
equation  (16)  gives  a  formula  for  *  in  terms  of  *  which  does  not 
contain  t  explicitly,  namely 


♦  = 


2(a2Q  +  2*)  (b20  +  2*) 

»  (»Z0  -  a20» 


% 

- 

< 

P  +  log 

b2o  +  2*  „ 

* 

- 

(17) 


Integration  of  equation  ( 17)  by  means  of  several  substitutions  yields 


where 


fe20  ‘  a  0  exp 

*  ~  t[  exp  (x)  -  lj  '  2 

b  n 

x  -£[1  -  exp  (kt/X)]  +  exp  (kt/X)  log  -j—  . 


(18) 


The  fields  of  radial  velocity  and  radial  and  circumferential 
strain  rate  can  now  be  written  down  from  equations  (4)  and  (5).  The 
stress  field  is  obtained  by  substituting  for  *  and  *  as  functions  of 
the  time  in  the  following  equations: 


4.  COMPARISON  OF  THE  VISCOUS  EFFECT  FOR  DIFFERENT 
VISCOSITY  COEFFICIENTS 

It  will  now  be  assumed  that  p  varying  with  the  time*  is  the 
internal  pressure  required  to  keep  the  expanding  tube  flowing  un¬ 
restrictedly  in  a  certain  manner.  In  order  to  simplify  the  problem 
the  condition  that  the  interior  boundary  of  the  tube  expands  with  a  un¬ 
iformly  increasing  speed  is  imposed.  The  pressure  p(t)  which  is 
required  to  produce  this  effect  is  found.  The  expansion  of  a  is  of 
the  form 

a  *  Sq  +  ct2,  (20) 


where  c  is  a  positive  constant, 
satisfy 

2  4  ? 

2*  =  c  t  +  2ayj  ct  , 

giving  *  =  2ct  ( a ^  +  ct2)  , 

••  2 
and  *  *  2c  (a^  +  3ct  ). 


The  function  *(t)  is  found  to 


(21) 


Omitting  details  of  the  calculation,  we  find  that  equation  (16)  now 
furnishes 


P(h)  =  £  log 


.  (-2» '  ,}1 

+  £ 

r  <- 

(3h,2+1)  los 

at  7  Z 

(1  +  h). 

+  F 

z  8 

1  + 


(«2o  -  »> 

d  +  h2) 


(1  -  «>2  I  1 


U-o>h 


(*20  +  hS  2h2)  (1+h2) 


(22) 


where  *q 
h 


Vv 

i 

t  (c/@q),  a  dimensionless  time. 
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and  P  (h)  =  p(t)/2k  ,  the  dimensionless  pressure.  Also  R  and  B 
are  two  dimensionless  parameters  analogous  to  the  Reynolds  number 
and  the  Bingham  number  and  defined  by 


For  the  present  purpose  it  will  be  more  convenient  to  utilize  less  signi¬ 


ficant  parameters  defined  by 


B  _  A- 


(23) 


The  pressure  can  then  be  written 


P(h) 


■  logL  t  S<s£u> 

2  L  (1  +h2)  J  p  2 


(«2o- 1) 

log  1  + - v- 2  + 

(1  +  h2) 


(1  -«0^k2 

(«  n  +  h4  +  2h^) 


-  B  (1-,<02)h 

P  /  4  2  2 

u  0  +h4  +  2h2)  (1+h2) 


Equation  (24)  gives  the  dimensionless  pressure  required  to  produce  an 
unrestricted  flow  of  the  tube  in  which  the  interior  boundary  has  a  uniform 
acceleration.  The  stress  field  is  then  given  by: 
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=  l  + 


log 


,  <  -“20> 

x  p 

P 

(3h2+l) 

(«20  +  h4  +  2h2) 

log 


i 

1 

H 

* 

.2,,  2,  ,.2  ,2  2  . 
h  (h  +1)  (p  q  -«  q) 

[  (o<2q  +  h4  +  2h2) 

(<4  +h4+  2h2)  (p2Q+h4+2h2) 

+  B,ji  ( 1  +h ^  (P20  +  *  0  +  2h4  +  4h2) 
U20  +  h4+  2h2)  (p20  +  h4  +  2h2) 


(25) 


and  the  velocity  field  u  is  found  to  satisfy 

u2  =  4  (  1  +  h2)2  a«ch2.  (26) 

- - -r-1 — y-  0 

(p2  +  h4  +  2h  ) 
where  pQ  =  r^a^. 

As  estimate  of  the  effect  of  viscosity  on  plastic  flow  of  the 
type  specified  above  was  obtained  from  equation  (24).  The  numeri¬ 
cal  values  used  in  the  investigation  were  as  follows  in  c.  g.  s.  units: 

^  =  5,*0=  2,  c  =  0.  1;  D  =  8.  5;  and  k  =  23.5  x  981  x  10^,  appropri¬ 
ate  for  a  thick  walled  brass  (Zn  30  percent*  Cu  70  percent)  tube.  The 
viscosity  coefficient  of  the  metal  in  these  units  would  be  of  order  10^. 

Two  values  X  =  5  x  10*®  and  X  =  10**  were  chosen*  and  the  results 
compared  with  that  for  X  =  0  corresponding  to  zero  viscosity.  With  the 
above  parameters  the  second  term  in  (24)*  the  inertia  term*  becomes 
comparable  with  the  other  terms  only  in  the  final  stage  of  yielding*  and 
for  all  practical  purposes  is  negligible.  The  variations  of  the  dimension¬ 
less  pressure  with  dimensionless  time* for  the  three  assigned  values 
of  the  viscosity*  are  shown  in  Fig.  1.  It  is  seen  that  the  viscosity  has 
a  considerable  effect  during  the  initial  stages  of  the  yield  process*  and  that 
in  order  to  maintain  the  same  yield  with  greater  viscosity*  the  initial  rate 
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of  loading  must  be  correspondingly  increased.  The  required  pressure 
attains  its  maximum  and  may  be  allowed  to  decrease  at  a  rate  more 
nearly  comparable  with  the  rate  of  unloading  in  the  nonviscous  case. 

In  each  case  the  required  pressure  approaches  zero  asymptotically 
with  increasing  time(  as  might  be  expected  in  the  absence  of  fracture* 
since  the  thickness  of  the  tube  is  steadily  decreasing. 


5.  COMPARISON  OF  PLASTIC,  VISCOUS  AND  INERTIA  EFFECTS 


It  is  also  of  interest  to  compare  the  effects  of  perfect  plasticity, 
viscosity,  and  inertia  on  the  expanding  tube.  For  this  purpose  an 
expansion 

a  =  (*q/4)  [  5  +  sin  {ut  -  rr/2]  (27) 


of  the  inner  radius  of  the  tube  is  imposed  during  the  time  interval 


t  =  0  to  t  =  ir/u  When  t  =  0,  a  =  aQ  and  da/dt  =  0;  when  t  =  (x/w), 
a  =  3a^/2  and  da/dt  =  0.  See  Fig.  2.  Hence  (27)  corresponds  to  an 
expansion  of  the  inner  radius  to  one  and  a  half  times  its  initial  value  in 
a  time  T  =  ir/  u),  starting  and  finishing  with  zero  velocity.  It  can  then 
easily  be  shown  that 


¥  = 


*  = 


*  = 


2 
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(-0 
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2 

a  0 

'  ~T~  • 
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/■ 
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a  o 

_  pnc  ?  j 

c— 0  ■ 
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(“•  'I 

16 

vUO  U  | 

32 

J 

(28) 

The  pressure  variation  needed  to  produce  expansion  (27)  can  be  written 
as 


Pjt)  =  Pp(t)  +  Pj(t)  +  Py(t)  ,  (29) 

2k 

where  P^  (t),  P^  (t)  and  P^(t)  are  the  contributions  due  to  perfect 
plasticity,  inertia,  and  viscosity  respectively,  and  given  by 
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15 


The  numerical  values  of  ag,  bg»  D  and  k  are  chosen  the  same 
as  previously.  The  value  of  X  is  now  fixed  at  5  x  10^.  By  varying 
«  the  expansion  can  be  made  to  take  place  in  any  desired  time.  Three 
values  of  w  were  chosen.  These  correspond  to  an  expansion  of  the  tube 
in  1  sec. »  10"  sec.  •  and  10  sec.  respectively.  The  pressure 

variation  with  time  required  to  produce  the  expansion  was  calculated 
in  each  case  from  equations  (28)  through  (30).  The  dimensionless  pres¬ 
sure  contributions  Pp,  Pj,  and  P^  were  plotted  against  t  in  each  case. 
When  T  =  1,  corresponding  to  a  slow  expansion  inertia  has  a  negligible 
effect  on  the  pressure  required,  but  the  perfectly  plastic  and  viscous 
contributions  to  this  pressure  are  of  the  same  order.  See  Fig.  3. 

When  T  =  10"^  the  required  pressure  is  considerably  increased  and 
the  perfectly  plastic  contribution  is  a  negligible  part,  but  the  inertia  ef- 

_7 

feet  begins  to  be  apparent.  See  Fig.  4.  For  T  *  10  the  inertia  term 
in  p(t)/2k  has  a  considerably  greater  effect  (see  Fig.  5). 


It  is  noted  that  Figs.  4,  5  show  that  a  negative  pressure  is  re¬ 
quired  near  the  end  of  the  period  to  produce  the  required  expansion.  By 
referring  to  equations  (14),  (15)  and  (16)  it  is  seen  that  p  can  be  replaced 

by  a  pressure  (p.  -  p  ),  where  p.  and  p  are  an  internal  pressure  and 
*  “  i  e 
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an  external  pressure  respectively.  Thus,  imposing  a  negative  value 
of  the  pressure  p  is  exactly  equivalent  to  an  application  of  external 
pressure  and  hence  is  physically  possible. 

The  order  of  magnitudes  involved  in  P  and  Pj  is  apparent 
from  equations  (28)  and  (30).  When  the  time  of  expansion  T  is  altered 
by  some  factor  K  ,  the  viscous  pressure  contribution  is  altered 
by  a  factor  K  and  the  inertia  contribution  by  a  factor  K  the 
perfectly  plastic  pressure  contribution  is,  of  course,  independent  of 
T.  For  the  type  of  problem  considered,  this  shows  an  interesting 
comparison  over  the  full  range  of  expansion  times  of  the  three  effects 
mentioned.  In  a  slow  expansion  the  flow  approximates  that  of  the  plastic 
quasi-static  theory  in  which  the  strain  rates  are  very  small  and  the 
inertia  effects  are  neglected.  As  the  speed  of  expansion  increases 
the  viscous  effect  becomes  important  and  then  dominant,  whereas 
the  inertia  effect  is  still  negligible.  For  even  faster  expansions 
the  effect  of  inertia  becomes  important  and  finally  predominates. 
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